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ABSTRACT 

The  acoustic  impedance  of  several  Helmholtz  resonators 
were  measured  using  two  different  techniques.  The  first 
technique  employed  the  conventional  Standing  Wave  Tube 
method.  The  sscond  was  the  Dual-Microphone  Transfer  Function 
method  of  J.Y.  Chung  and  D.  A.  Blasec .  The  calculahility  of 
the  frequer.cy  dependent  Helmholtz  resonator  impedance 
allowed  the  methods  to  be  compared  for  both  absolute  accu- 
racy and  relative  precision.  The  average  relative  precision 
for  the  SWT  technique  was  4.2%  while  the  average  for  the 
dual-microphone    technique   was    19.4%. 
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I-    INTRODUCTION 

Acoustic  impedance  tubes  have  been  used  for  many  years 
to  obtain  complex  acoustic  impedances,  reflection  coeffi- 
cients and  phase  information  for  nan/  different  types  of 
materials.  The  techniques  used  with  the  impedance  tubes  have 
been  many  and  varied.  This  study  involves  the  comparison  of 
the  dual-microphone  transfer  function  technique  used  by 
Chung  and  Blaser  [Ref.  1]  and  the  standing  wave  tube  (SWT) 
technique  to  determine  their  absolute  accuracy  and  relative 
precision.  This  work  was  undertaken  as  a  prelude  to  the 
extension  of  computer  controlled  impedance  measurement  tech- 
niques  to   water. 

A.       IMPEDANCE    MEASUREMENT    TECHNIQUES 

Many  different  impedance  measurement  techniques  have 
been  used  over  the  years.  Beranek  [Ref.  2]  divides  the 
different  techniques  into  three  general  categories,  Surface 
Methods,  Transmission  Line  Methods  and  Comparison  Methods. 
This  study  looked  closely  at  only  two  of  the  techniques, 
however,  seme  of  the  other  methods  are  worth  mentioning  as 
background. 

1  •      Surface    Methods 

These  methods  involve  the  measurement  of  the  sound 
pressure  at  a  point  on  the  surface  of  interest  and  the 
particle  or  volume  velocity  at  the  surface.  Cook  [Ref.  3] 
did  this  using  a  short  tube  of  a  given  radius  and  length  and 
a  clamped-edge  diaphram.  He  was  then  able  to  calculate  the 
volume  velocity  and  the  average  pressure  over  the  surface  of 
the    sample.      Beranek   [Ref.    4]   presented      a    useful    theory   for 


predicting  the  normal  acoustic  impedance  of  porous  materials 
in  terms  of  their  flew  resistance,  porosity  and  effective 
dynamic    mass   of    the    enclosed   gas. 

2-      Transmission    Line   Methods 

By  f ar ,  the  most  popular  methods  of  impedance  meas- 
urements fall  into  this  category.  The  standing  wave  tub?  and 
daai-microphone  techniques  in  this  study  are  included.  They 
will  be  examined  in  detail  later.  Another  technique  involves 
the  measurement  of  the  pressure-frequency  resonance  curve 
and  determining  the  curve  width  and  resonant  frequency. 
Harris  [Ref.  5]  [Ref.  6]  discusses  two  of  these  methods 
which   are   closely   related   to   each   other. 

Finally,  one  additional  method  currently  in  use  by 
the  U.S. Navy  for  making  acoustical  measurements  deserves 
mentioning.  [Ref.  7]  This  system  is  digital  and  is  used  for 
underwater  acoustical  measurements  from  1  Hz.  to  approxi- 
mately 2  MHz.  This  digital  measuring  system  employes  a 
buffered  video  analog-to-digital  converter  and  was  primarily 
designed  to  perform  acoustical  measurements  using  the  Prony 
Method.  [Ref.  8]  [Ref.  9]  The  system  uses  a  stable  signal 
and  sampling  source  and  12-bit,  5  MHz.,  analog-to-digital 
converter  with  a  UK  buffer  memory.  A  PDP-11/23  serves  as  a 
controller  as  well  as  a  signal  processor  for  averaging  and 
amplitude   and   phase    estimation. 

3.      Comparison    Methods 

These  methods  mainly  involve  the  use  of  acoustical 
bridges  or  reliable  acoustical  standards.  Steward  first 
introduced  [Ref.  10]  the  simplest  and  most  rapid  to  use. 
Later  refinements  followed  [Ref.  11]  but  suitable  variable 
standards  of  acoustical  resistance  and  reactance  were  diffi- 
cult to  obtain.  The  use  cf  reliable  acoustic  standards, 
rigid    wall    and      an    eighth-wavelength    tube,       was      reported   by 
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Flanders  [Ref.     12]    in    1932.       However,      mathematical   calcula- 
tions  were   difficult. 

In  'his  study,  the  transmission  line  method  and  the 
comparison  method  are  combined  by  introducing  a  new  "stan- 
dard", the  Helmholtz  resonator.  This  standard  has  the 
advantage  of  having  an  easily  calculable  reactance  which  is 
absolute  to  within  one  adjustable  parameter.  This  parameter 
is  fixed  by  the  Kelmhoitz  resonance  frequency  which  is 
easily  measured  to  a  high  degree  of  precision  and  absolute 
accuracy. 

B.       HISTORY    OF    DUAL- MICROPHONE    MEASURING    TECHNIQUES 

Clapp  and  Firestone  [Ref.  13]  first  mentioned  the  use  of 
two  microphones  in  an  article  published  in  1941.  In  this 
article  it  was  recognized  that  their  equipment  could  be  used 
to  measure  the  absorption  coefficient  of  materials.  Later, 
Shultz  [Ref.  14]  published  an  article  on  the  acoustic  watt- 
meter and  its  application  for  measuring  absorption  coeffi- 
cients. During  this  time  the  use  Df  the  acoustic  wattmeter 
for  impedance  measurements  was  little  used  since  the 
standing  wave  tube  (SWT)  technique  was  the  most  favored 
method. 

In  1977  Seybert  and  Ross  [Ref.  15]  published  an  article 
which  discussed  an  experimental  technique  using  two  micro- 
phones and  random  noise  excitation.  This  technique  used 
digital  methods  and  the  auto-  and  cross- spectra  to  determine 
the  normal  complex  reflection  coefficient  cf  a  sample  termi- 
nating a  tube.  Finally,  in  1980,  3hung  and  Blassr  [Ref.  1] 
published  their  article  discussing  the  transfer  function 
technique  for  measuring  in-duct  acoustic  properties.  This 
technique  also  used  digital  methods  and  the  auto-  and 
cross-spectra  to  determine  the  normal  complex  reflection 
coefficient.      However,      the    mathematics      did    nor    involve   the 
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matrix-inversion    as  =  d      by  Seybart   and      Ross    and   as      a    result 
the    ccmput stior.    time    was   significantly   reduced. 

A  two  microphone  technique  utilizing  pure  tone  rather 
than  random  excitation  and  analog  rather  than  digital 
processing    is   described   by    Elliot.      [Ref.    16] 

C.       PDSPCS2    OF    STUDY 

The  above  methods  have  been  used  to  measure  numerous 
materials  of  differing  acoustic  properties.  However,  the 
properties  of  these  materials  were  not  known  or  calculable 
to  any  reasonable  accuracy.  The  aim  of  this  study  was  to 
determine  the  absolute  accuracy  of  the  dual-microphone  tech- 
nique and  its  relative  precision  compared  to  the  SWT  tech- 
nique as  a  prelude  to  the  development  of  an  impedance 
measurement  apparatus  for  use  in  water.  A  Helkholtz  reso- 
nator   was   used    to   accomplish   this   comparison. 
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II-  COMPARISON  QF  MEASUREMENT  TECHNIQUES 

Previous  studies  of  impedance  measurement  techniques 
involved  the  use  of  materials  for  which  the  impedance  was 
not  accurately  known.  This  study  deals  with  impedance  meas- 
urements using  a  Helmhoitz  resonator  as  a  calculable  stan- 
dard to  determine  the  absolute  accuracy  and  relative 
precision  of  the  two  techniques.  A  Helmhoitz  resonator  was 
chosen  bacause  its  various  acoustical  parameters  could  be 
independently  measured  or  calculated  and  then  compared  to 
the  experimental  results. 

A.   THEORETICAL  BACKGROUND 

The  term  "acoustic  impedance"  oan  have  different  mean- 
ings for  different  situations.  It  is  necessary  to  define 
exactly  the  different  impedances  involved  in  this  study. 
This  is  necessary  to  prevent  confusion  later  when  the  rela- 
tionship between  the  different  impedances  is  used  to  deter- 
mine the  calculated  acoustical  properties  of  a  Helmhoitz 
resonator.  Normally  this  is  not  necessary  since  the  specific 
acoustic  impedance  and  the  acoustic  impedance  differ  by  a 
constant  with  the  dimensions  of  area. 

1  •   Specif  j.c  Acoustic  I  mpe dance 

The  specific  acoustic  impedance  is  defined  as  the 
ratio  of  acoustic  pressure  in  a  medium  to  the  associated 
Eulerian  particle  speed.   [  Ref .  17] 

Z5  =  p/u  (2.1) 

The    MKS    unit   of    specific   acoustic   impedance    is   Pa-s/m. 
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The  specific  acoustic  impedance  of  a  medium  is  a 
real  quantity  for  plana  waves.  H:>;/ever,  this  net  true  fer 
the  standing  waves  generated  in  an  imps  dance  tube.  ?c:  this 
casa,      the   specific      acous~ic    impedance    will    be      complex   and 

can    expressed   in   the    form 


Z.s    =    r   +    jx 


(2.2) 


where    "r"    is  the   specific  acoustic   resistance    and    "x"    is   the 
specific  acoustic   reactance    for   the    medium    concerned. 

2.      Acoustic   Imogdance 

The  acoustic  impedance  of  a  medium  acting  en  a 
surface  of  area  "A"  is  the  ratio  of  the  acoustic  pressure  at 
the    surface   to    the    vclume  velocity   at    the   surface    [Ref.    17] 


Z*.    =    P/U 


(2.3) 


where  U  =  uA  if  u  is  normal  to  the  surface  of  area  "A".  The 
HKS  unit  of  acoustic  impedance  is  Pa-s/ra3,  which  is  also 
termed  an  acoustic  ohm.  The  acoustic  impedance  and  the 
specific  acoustic  impedance  are  related  at  the  surface  of 
area  "A"  by 


Z-  =  Z./A 


(2.4) 


3  .       Mechanical    I  mpedar.ee 

The  mechanical  impedance  is  the  ratio  of  the  driving 
force  to  the  resultant  particle  speed  at  the  point  where  the 
force    is   applied  [Ref.    17] 


Z„   =    F/u 


(2.5) 


The    MKS    unit  of    mechanical    impedance    is    N-s/m,    also    referred 
to    as    a    mechanical    ohm. 
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y •      Helmholtz  R^scn a tor 

A  Helmholtz  resonator  can  be  treated  as  a  lumped 
pax  am  ST.  3  ;r  system*  Ine  risen  ax  or  consxsts  or  a  vcxum  = 
with  rigid  walls  and  a  neck  of  length  "L"  and  cross- 
sectional  area  "A".  The  acoustic  pressure  in  the  volume 
provides  the  lumped  stiffness  "s"  if  the  wavelength  is  much 
greater  than  the  cube  root  cf  the  volume.  The  fluid  in  the 
neck  provides  the  lumped  mass  as  long  as  the  wavelength  is 
much  greater  than  the  length  of  the  neck.  Finally.  the 
resistance  is  provided  by  the  opening  radiating  sound,  if 
the  wavelength  is  much  greater  than  the  sguare  root  of  the 
cross-sectional  area  of  the  neck.  An  additional  resistance 
is  introduced  by  viscous  losses  in  the  neck  and  irreversible 
thermal  conduction  to  the  volume  boundries  if  the  compres- 
sion of  the  gas  is  adiabatic  and  the  walls  remain  isoth- 
ermal. However,  since  these  resistance  sources  are  not 
easily      calculated         from      "First     Principles",  they      are 

neglected.  The  radiation  resistance  "R  "  can  be  calculated, 
however,  its  value  is  dependent  upon  the  impedance  of  the 
medium  which  changes.  To  get  the  viscous  loss  resistance 
term  the  flow  field  in  the  neck  is  required.  Therefore  these 
resistance    terms  are   cf   little   use   in    this    calculation. 

The    fluid   in    the   neck    has   an    affective   mass    given   by 
[Ref.    17] 


m    =  p  AL ' 


(2.6) 


where  "L'  "  is  the  effective  length  of  the  neck  and  is  a 
function  of  actual  length  "L"  and  the  neck  radius  "a"  and 
the  flew  field  in  the  vicinity  of  the  oriface.  The  effective 
mass  is  greater  than  the  actual  mass  in  the  neck  due  the 
radiation-mass  loading  at  the  ends  generated  by  this 
extended  flow  field.  The  stiffness  of  the  resonator  is 
given   by  [Ref.     17] 


s    =     (/3CZ)  (A*  )/V 


(2.7) 
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By  driving  the  Helmholtz  resonator  with  a  sound  wave 
of  pressure  amplitude  "?",  the  following  equation  for  th<= 
inward  displacement  "x"  of  the  fluid  in  the 
[Ref.  17] 


nark     T"  -•  ^"1  1  t  - 


m^~*    +    R^"^    +    sx    =    APSW 


(2.8) 


This  differental  equation  is  the  same  as  that  fcr  a  driven 
oscillator  and  the  irechanical  impedance  of  the  Helmholtz 
resonator  can   be   expressed    as   [Ref.    17] 

Z„    =    Rr    +    J  (om   -    S/fe)  (2.9) 

where   the    mechanical    reactance    "X^"    is   given    by 

X„,    =  torn    -  s/co  (2.  10) 

Since  the  resonance  occurs  when  \m  -  0,  the  resonant 
frequency   nui  "    follows 

u>*  =   s  /m  (2.11) 

Using   equations     (2.6)    and    (2.7) 


The    final   result   obtained   from   cancelling   terms    is 


to*  =    c 


/^F 


(2. 12) 


It  should  be  noted  that  nowhere  in  the  derivarion  of 
equation  (2.12)  is  any  restriction  placed  on  the  physical 
arrangement  of  the  Helmholtz  resonator.  Figure  2.1  shows  the 
arrangement      used      for  this      study.         As      can    be      seen      from 
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equation  (2.12)  and  the  figure,  the  resonant  frequency  can 
be  changed  by  changing  -ha  area  of  "neck"  (diameter  D,  ),  -he 
volume  of  the  resonator  (change  depth  h)  or  the  length  of 
the    "reck"  (thickness    t) 

Recalling  equation     (2. 10) 


Xm    =  u)m    -   s/£j 


(2. 10) 


and    by      factoring  out      "a"    and      using    equation      (2.11)  f      the 
following   result   is    obtained 

x„,    =   m(u>  -  rffo) 

Since    io  =   2rrf  the   final   form   is   obtained 
X^    =    2/rm(f   -    C/f) 

where    f   is   the    frequency   in    hertz.    For   convenience    we   define 
X«    =    f    -    f|/f  (2.  13) 


so    that 


X„    =    2^mX«, 


(2.14) 


I*      is      now      possible      to      derive     the     relationship 
between    the    specific    acoustic   reactance    "Xy  "    and    "X^".    Given 

a  resonator   as    shown    in    Figure   2.2,       divided    into    sections    1 
and    2    as   indicated    and   recalling    equation     (2.1) 


ZS    =    Pi  /»/ 


(2.15) 


where   p;    and     u,    represent    the    pressure      and    particle    speed, 
respectively,      in   section    1 .      This      is   the    specific    acoustic 
impedance   of  the   medium   in    section    1. 
Next  asing    equation     (2.5) 


lm   =    F*/Uz 


(2. 16) 


Converting    F,    to   a    pressure    at   the   boundary   results   in 
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Z„  =    -&A-  (2.17) 


where  pr  and  u,  represent  the  pressure  and  partial-  steed, 
respectively,  in  section  2.  Az  is  the  cross-sect ic nal  area 
of  the  resonator  neck.  Since  there  is  no  physical  boundary 
between  sections  1  and  2,  the  pressures  on  both  sides  must 
be    equal. 

P,     =    Fz    =  P 

Also,  the  volume  velocities  must  be  equal  across  the 
boundry.    Therefore 

uz     =    -^'-u,  (2.  18) 

Substitutinq    equation       (2.13)       into      equation    (2.17) 
and    reqrcupinq    the    appropriate  terms    results    in 

Zm     =       (-£-)    ("$-)  (2.   19) 

Finally,  it  can  be  seen  that  the  quantity  (p/u,)  is  just  the 
specific   acoustic   impedance    and  the   result    is 


~2l~    =   ~%c  (2.20) 


Since  the  impedances  are  composed  of  the  resistive  and  reac- 
tive parts,  the  reactive  portions  vary  in  the  same  manner. 
Therefore, 


•£-    =    A'*  (2.21) 
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It    car.    be   seen  from    equation     (2.14)    -chat 


-H   =    2/7m~&- 


x:  x 


»r* 


'V\ 


It  should  be  noted  that  the  two  techniques  in  this 
study  provide  the  values  for  specific  acoustic  impedar.cs 
normalized  to  the  characteristic  impedance  "  />  c" .  Therefore, 
the  preceding  equation  becomes 


and,  finally,   substituting  for  "m"  from  equation  (2.6)   and 
for  "X^/X*,  "  from  equation  (2.21) 


'xf  =  ZILc  '#?'—  (2*22) 


As  can  be  seen,  this  equation  does  not  depend  on 
frequency  and  therefore  remains  constant  as  frequency 
changes.  For  convenience  this  constant  is  called  "Ky^"  with 
the  subscript  signifying  the  theoretical  value  for  the  ratio 
of    "XyZ/scX^,"   for   a    given   Helmholtz   resonator. 

Km  =    ^-—jr-^ (2.23) 


L 


"Ky.j-  "  has  been  reduced  to  an  expressicn  which 
contains  only  directly  measurable  terms  with  the  exception 
of  "L*"r  the  effective  length  of  the  neck.  However,  by  using 
equation  (2.12)  "L ' "  can  be  determined  from  parameters  which 
are  measurable  by  separate  means.  Later  in  this  chapter, 
"L* "  will  be  looked  at  more   closely  for  comparison  with  the 
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calculated    "L"' 
configurations . 


for 


ether 


Helmholtz 


resonato: 


5. 


STANDING  WAVE  TUEE  TECHNIQUE  (SWT) 


The  standing  wave  tube  technique  provides  a  convenient 
laboratory  method  for  determining  the  absorption  coefficient 
and  specific  acoustic  impedance  of  materials.  The  method, 
however,  is  somewhat  tedious  when  a  large  number  of  samples 
is  measured  or  measurements  of  a  single  sample  are  required 
at.  a  great  number  of  frequencies. 

1 .   Theory 

The  sample  to  be  measured  is  placed  at  one  end  of 
the  standing  wave  tube.  At  the  other  is  mounted  a  loud- 
speaker. A  microphone  probe  is  inserted  through  the  loud- 
speaker end  and  is  attached  to  a  microphone  car.  The  car 
moves  such  that  the  microphone  probe  is  moved  along  the  axis 
of   the   tube. 

A  pure  tone  signal  is  reflected  from  the  sample.  The 
reflected  wave  interferes  with  the  incident  wave  and  a 
standing  wave  is  established.  The  microphone  probe  and  car 
arrangement  allows  fcr  the  location  and  acoustic  pressure  of 
the  nodes  and  anti-nodes  to  be  measured.  From  the  measure- 
ment of  these  amplitude  and  phase  characteristics,  calcula- 
tion of  absorption  coefficient  and  specific  acoustic 
impedance   of  a    sample  is   possible. 

From  reference  6,  the  following  are  the  equations 
needed  for  the  determination  of  the  desired  acoustic 
parameters. 


X 


(2.24) 


where 
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-; 


^  =    coth 


(2.25) 


and 


1 f~ 


(2. 26) 


.CC^l-iC 


n       J-^i 


np  m  aid  "P  "  are  -v 
standing  wave  tub*  and  "d"  is  the  distance  to  the  firs- 
pressure  minimum  (node)  from  the  sample.  (NOTE:  The  micro- 
phone output  voltage  is  the  actual  measurement  made  since 
only    the   ratio    of   the    pressures    is   required) 

2.      Equipment   Set-up 

Figure  2.3  is  a  simplified  representation  of  the 
equipment  used  in  the  standing  wave  tube  technique.  It 
includes  the  Bruel  and  Kjaer  Standing  Wave  Apparatus  Type 
4002,  a  HP  302A  Wave  Analyzer  and  an  amplifier  and  frequency 
counter. 

a.      3ruel      and        Kjaer      Type         4002      Standing        Wave 
Apparatus 

The  apparatus  is  designed  for  easy  and  quick 
evaluation  of  acoustic  materials  by  the  standing  wave 
method.  A  major  advantage  is  that  relatively  small  circular 
samples    (10   cm.    in    diameter)    are   all    that    are   required. 

The  tube  used  had  an  interior  diameter  of  10 
cm.,  which  covered  a  frequency  range  of  90  to  1800  hz.  The 
lower  limit  was  set  by  the  length  of  the  tube  and  the  upper 
limit  determined  by  the  cut-off  frequency  for  the  lewest 
non-plane      wavemode      of      the   tube.  [Ref.    18]      Two      sample 

holders  were  also  used  with  fixed  depths  of  1  in.  and  2  in., 
respectively.  The  bases  of  the  holders  were  very  thick  to 
present    as    rigid   a    termination    as    possible. 


21 


The  microphone  in  the  microphone  car  is  a  small 
dynamic  microphone  -of  adequate  frequency  range  for  the 
standing  wave  tube  method.  It  is  isolated  by  elastic 
supports  in  the  car  to  reduce  the  effaces  of  external 
airborne  noise    and    impact   sound   or    vibration. 

b.      Hewlett-Packard    HP    302A    Wave    Analyzer 

The  Hewlett-Packard  medal  302A  wave  analyzer  is 
a  tunable  voltmeter  cf  high  sensitivity  over  the  frequency 
range  of  20  Hz.  to  5  0  kHz.  The  wave  analyzer  in  the  3F0  mode 
of  operation  becomes  an  audio  signal  generator  combined  with 
a  tuned  voltmeter.  Further,  since  the  tuned  voltmeter  has  a 
very  narrow  pass  band  (6  Hz),  noise  and  harmonics  of  the 
signal   have   very   little   effect   on   the    voltmeter   reading. 

3  .      Experimental   Proced  ures 

The  standing  wave  tube  technique  was  used  in  several 
different  experiments.  It  was  used  with  a  rigid  end  as  the 
sample  to  measure  the  "acoustic"  end  of  the  microphone  probe 
tube  and  the  sound  speed.  It  was  also  used  to  determine  the 
resonant  frequency  and  reactive  impedance  of  three  Helmholtz 
resonators.  Table  I  shows  the  physical  dimensions  of  each  of 
the    resonators. 

a.       Data   Collection 

In  all  of  the  experiments  the  data  was  collected 
in  the  same  manner.  However,  in  the  first  experiment  (riqid 
end)  the  only  data  collected  was  the  location  of  the  first 
two  nodes  and  the  frequency  of  the  signal.  Table  II  shows 
this    data. 

For  the  experiments  involving  the  Helmholtz 
resonators  the  values  of  the  maximum  and  minimum  microphone 
output  voltages  were  obtained  in  addition  to  the  frequency 
and    node   locations.    This   data    was   placed    in   Tables    III-V 
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TABLE    I 
H^lmholtz    Resonator   Physical   Dimensions 


Besona 

tor 

Volmue 

N9 

ck    Diam 

eter 

Neck   Length 

no. 

cm  . 

cm. 

cm. 

1 

195.52 

2.52 

1.27 

2 

381  .01 

2.52 

1.27 

3 

195.52 

3.46 

1  .00 

b.       Data    Reduction 

In  the  first  experiment,  the  measured  distance 
between  the  two  nodes  was  used  to  calculate  the  wavelength. 
From  this,  the  quarter- wavelength  was  determined  and 
compared  to  the  distance  from  the  sample  to  the  first  node 
as  measured  by  the  microphone  car.  Since  the  "sample"  was  a 
rigid  end  located  5.00  cm.  from  the  measured  end  of  the 
tube,  the  two  distances  should  have  differed  by  5  cm. 
Additionally,  the  wavelength  and  the  frequency  were  used  to 
calculate  the  sound  speed  in  the  tube  for  each  frequency 
measured.  Thsse  results  are-  shown  in  Table  II  The  average 
value  of  "  d"  suggests  that  the  "acoustic  center"  of  the 
probe  tube  is  1.7  +  .5  mm.  ahead  of  the  physical  end  of  the 
probe    tube   over    the    frequency    range   examined. 

For  the  experiments  involving  the  Helmholtz 
resonators,  the  reduction  of  the  data  was  somewhat 
different.  The  wavelength  was  determined  by  one  of  two 
methods      depending      upon      whether   two      node      locations      were 


23 


TABLE    II 
Stan  ding    Have  I 


.  U  &        fcj   \J  < 


Freg.        First   Min.        Sec.    flin.         (*£tL)  ad               c 

Hz.                  cm.                          cm.                  cm.  cm.  cra./s. 

299                23.90                      81.74              28.93  5.03  34588 

400                  16.80                       59.95               21.58  4.78  34520 

600                   9.58                      38.38               14.40  4.82  34560 

800                    5.96                       27.53               10.78  4.79  34512 

1001                    3.84                       21.03                 3.60  4.76  34414 

1197                    2.19                       16.91                 7.35  5.16  35240 

1402                     1.40                        13.64                  6.13  4.73  34321 

1599                    0.54                       11.36                 5.40  4.86  34602 

average  4.37  34595 

deviation  .150              277 

per    cent  03.08  0.80 


i 


measurable.  In  those  cases  where  two  or  more  node  locations 
were  measurable,  the  distance  between  them  was  used  to 
calculate  the  wavelength.  At  the  lower  frequencies  where 
only  one  node  was  measurable,  the  wavelength  was  determined 
from  the  frequency  and  the  average  sound  speed  found  during 
the  experiment.  It  should  be  noted  that  an  attempt  was  made 
to  use  the  location  cf  the  measurable  anti-node  to  determine 
wavelength.  However,  this  measurement  was  too  inaccurate  due 
to  the  nature  of  the  waveform  being  measured.  At  the 
maximum,  the  values  of  amplitude  are  changing  very  slowly 
(second  order)  so  that  the  measurement  of  the  exact  location 
of  the  maximum  was  very  difficult.  On  the  other  hand,  at  the 
nodes,  the  values  were  changing  more  rapidly  (first  order) 
and   the    location   could  be   easily   and    precisely  determined. 
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TABLE    III 
Standing    Wave    Tuba   Measured   Data   -    Resonator    #  1 

Freq.  First  Min.  Sec.      gin.  Third  Bin.         Maximum 

Hz.  cm.  uv.  cm.         uv.  cm.  uv.  mv. 

250  31.75  65  2.20 

350  20.8  5  5  0  69.70         6  0  1.30 

400  16.65  260  59.69      260  4.60 

450  12.30  140  50.40       130  39.20  175  1.80 

500  6.70  450  41.20      460  76.30  520  2.80 

550  1.70  240  30.35       290  61.65  190  1.60 

600  .55  200  22.50       180  51.30  185  1.78 

650  .10  200  18.42       140  45.10  165  2.60 

700  15.90         34  40.64  31  1.00 

750  14.20         90  37.15  100  5.20 

1000  9.75         30  27.15  20  3.40 

1250  7.60            5  21.40  5  1.95 

1500  6.75         25  17.80  5  2.50 


TABLE    IV 
Standing    Wave    Tube   Measured   Data   -   Resonator    #2 

ZL§SL-         First  Min.  Sec.  Mi n .  Third  Min.         Maximum 

Hz.             cm.  uv.  cm.  uv.  cm.  uv.              mv. 

250  27.88  1750  22.1 

300  19.69  335  76.72  330                                              2.90 

350  10.12  1610  59.32  1870                                                 8.45 

375            3.65  1280  49.63  1510                                                 6.70 

400            1.59  800  40.94  915  84.65  355               4.35 

450              .43  480  29.90  615  68.58  565               8.00 

500            0.00  164  23.96  75  58.64  128               2.88 

750  13.53  24  36.35  52  59.50  38               2.91 

1000            9.75  13  26.96  15  44.33  39               2.43 

1250            7.44  24  21.71  43  35.12  88               2.00 

1500            6.32  30  17.93  10  29.37  42               4.42 
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TABLE    V 

St 

an ding    Wave    Tub 

3  Measured   Da 

ta   ~   3 e so sat o: 

-    fi  3 

Fr  ec[. 

First 

Hin- 

Sec . 

Bin. 

Third 

Hiii. 

Maximum 

HZ. 

cm. 

mv . 

era . 

mv. 

cm. 

mv. 

mv. 

200 

41.61 

3.20 

51  .5 

250 

31.74 

29.  5 

5  97 

300 

26.20 

2.90 

83.  95 

3.  15 

55.5 

350 

21.51 

3.90 

70.  04 

3.30 

65.5 

400 

18.42 

10.  0 

61  .  49 

9.40 

178 

450 

16.17 

1.  52 

54.  44 

1.70 

26 

500 

14.26 

4.05 

48.  19 

3.40 

82.21 

4.  05 

61.5 

550 

1  1.10 

5.40 

42.  92 

5.30 

73.68 

5.60 

68.5 

600 

9.18 

2.05 

37.66 

1.98 

66.69 

1.70 

24.0 

650 

6.43 

6.  05 

32.  79 

6.  15 

60.39 

6.  10 

46.  2 

700 

4.28 

4.00 

28.  62 

3.95 

52.  94 

3.80 

26.3 

750 

2.29 

2.90 

24.  54 

2.55 

47.49 

2.  64 

15.4 

800 

19.07 

3.95 

40.  41 

4.25 

62.  37 

4.  15 

27.3 

850 

16.52 

1.10 

37.  32 

1.  12 

57.42 

1.  10 

10.  3 

900 

14.41 

3.10 

33.  30 

2.90 

52.80 

3.20 

47.  5 

1000 

1  1.67 

.80 

29.  00 

.90 

46.  18 

.76 

22.5 

1250 

8.22 

.60 

21  .  86 

.45 

36.03 

.52 

24.6 

1500 

6.53 

.65 

17.  96 

1.40 

29.49 

1.  13 

60.5 

After  the  wavelength  was  determined,  it  was  used 
to  scale  the  distance  to  the  first  node  (Table  VI)  .  These 
values  were  plotted  versus  the  frequency  (Figures  2.4-2.6). 
The    first   two   points    above    resonance    were   also  plotted    as 

d/A     -    .5  (2.27) 

A  linear  regression  was  then  used  with  the  two  data  points 
both  above  and  below  resonance  to  determine  the  resonant 
frequency  "f  ".  The  resonant  frequency  for  each  resonator  is 
also    shown    on   Figures   2.4-2.6. 

Once  the    resonant    frequency      was   obtained,      "X,J," 
was   calculated    for    each   frequency   using    equation     (2.13).      "X, 
//jc"      was      then      calculated      using      equations       (2.24-2.26). 
Finally,      the   ratio    of   "Xy/zsc"    to    "XJJ"     (K^p)       was    calculated 
for    comparison    to  the  theoretical   value    (K^)     determined   from 
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TABLE    71 

Standing  Wave   Tube  i/^j     Data 

Freer.  Pes on at or    #1  R  e senator    #2        Resonator    #3 

hz.  a /a  a/*  a/A 

200  .2416 

250  .2437  .  2015  .2303 

300  .220  1  . 1726  .2268 

350  .21  34  .  1028  .2214 

375  . 0397 

400  .1934  .  4683(-.0317)        .2138 

450  .1585  .  3865  (-.11  35)        .2113 

500  .0954  .  3454  .2099 

550  .4848  (-.0152)  .1774 

600  .3906(-.1G94  .1596 

650  .3452  .1192 

700  .3213  .0880 

750  .3094  .2943  .0507 

800  .4404(-.0596) 

850  .4039(-.0961 

900  .3754 

1000  .2802  .2820  .3382 

1250  .2754  .2688  .2950 

1500  .2706  .2751  .2844 


equation    (2.23).       Tables    VII-IX    list    the    experimental    values 
for    each   Helmholtz    resonator. 

c.      Error   Analysis 

Before  an  analysis  of  the  errors  could  be 
attemptea,  it  was  necessary  to  identify  all  of  the  sources 
of  the  errors.  The  errors  were  a  result  of  the  measurements 
of  the  voltage  (P,,,^  ar.a  P^/0)  ,  frequency  ana  aistances  to  the 
nodes. 

Recalling  equation     (2.25) 


V,=   coth"'    (-§»*«---)  (2.25) 
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TABLE    711 
Standing    tfave   Tube   Reduced   Data   -    Resonator    #1 


Frequency 

c 

Hz. 

cm.  /s 

250 

32050* 

350 

34  195 

400 

34  43  2 

450 

34  60  5 

500 

34  80  0 

550 

34430 

60  0 

34560 

650 

34  6  34 

700 

3463  6 

750 

34  42  5 

1000 

34800 

1250 

3450  0 

1500 

34  6  50 

average 

34560 

deviation 

174 

per  cant 

.50 

f„=    542.7 

Hz. 

Is^oc 

X  • 

Q^JCP 

Pa-s/m 

S  l                   X 

10-3 

-7     3Q 

-928. 1 

7.  °6 

-4.  15 

-491.5 

8.44 

-2.63 

-336.3 

7.82 

-1.  52 

-204.5 

7.43 

-.658 

-39.05 

7.  39 

.094 

14.50 

6.48* 

.807 

109.  1 

7.40 

1.45 

196.9 

7.36 

2.07 

279.3 

7.41 

2.55 

357.3 

7.  14 

5.  21 

705.5 

7.38 

6.22 

1014 

6.  13 

7.70 

1304 

5.90 

average 

7.57 

deviation 

.39 

per   cent 

5.  15 

*    Values  not    used  for   the   calculation   of   the    average, 


it      can    be     seen  that      for    the      values   of      "PmM"      and    nVnt^  " 

involved,      SWH        is    small      and    therefore,        an   error      in   the 

p 
standing   wave  ratio     (SWR   =      mjf~Sm   I    has    1^'ttl9  effect,   on   the 

overall   calculations. 

Looking    at   equation    (2.24) 


and    taking    into    account   that   "y"    is    small,      equation    (2.24) 
then    reduces  to 
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TABLE   VIII 

i 

Standing 

Wave    Tube 

Reduced   Data   -    Resonator    *2 

f,=    387.8 

Hz. 

Freguer. 

SI 

c 

X,^c 

£**r 

Hz. 

cm. /s. 

Pa-s/m 

S_1                X 

10-3 

200 

34392* 

-3.86 

-551.9 

6.99 

250 

34430* 

-2.43 

-351 . b 

6.91 

300 

34  54  2 

-1.46 

-201.3 

7.25 

350 

344  4  7 

-.563 

-79.69 

7.06 

375 

33563* 

-.070 

-26.04 

2.69* 

400 

34296 

.258 

24.03 

10. 74* 

450 

3351  6 

.745 

115.8 

6.43 

500 

34  62  0 

1.55 

199.2 

7.78 

550 

34557 

2.08 

276.6 

7.52 

600 

34  39  8 

2.35 

34  9.4 

6.73 

650 

34437 

2.  30 

418.6 

6.69 

700 

3450  3 

3.23 

485.2 

6.66 

750 

34508 

3.72 

549.5 

6.77 

1000 

34560 

5.49 

849.6 

6.46 

1250 

349  13 

8.51 

1130 

7.53* 

1500 

34  4  70 

7.  11 

1400 

5.08* 

average 

34  44  4 

avera  ga 

6.  94 

deviat 

ion 

314 

deviation 

.  41 

per   c 

ent 

.91 

per   cent 

5.9 

♦Value 

s   not 

used    fcr 

the   calcula-1 

:ion    of   the 

average. 
....   _j 

or 


-yt-  *  cot  ^ 


away    from   resonance. 

On  the  basis  of  this  result,  it  was  determined 
that  the  best  data  was  obtained  in  the  regions  just  away 
from  resonance  for  the  following  reasons.  First,  since  tar- 
goes  tc  zero  at  resonance,  "  <£."  terms  begin  to  dominate  and 
the      errors   in       "y, "      become   large      when      compared      to    n  <f>^* . 
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TABLE    IX 
Standing    Wave   Tube  Reduced   Data   -   Resonator    #3 


f„=    773.0 

Hz. 

Freauency_ 

c 

Xs4°c. 

X  ' 

Zejip 

Hz. 

cm. /s. 

Pa-s/m 

s-i 

X 

10-3 

200 

-7.86 

-2788 

2.81 

250 

-6.93 

-2140 

3.24 

300 

34650 

-5.72 

-1692 

3.  38 

350 

339  9  9 

-5.03 

-1357 

3.71 

400 

34  45  6 

-4.09 

-1094 

3.74 

450 

34  44  3 

-3.78 

-    878 

4.  31 

500 

33975 

-3.  63 

-    695 

5.22* 

550 

34  41  9 

-2.07 

-    536 

3.  86 

600 

34  50  6 

-1.50 

-    396 

3.79 

650 

35074 

-.870 

-    269 

3.  23 

700 

34062 

-.600 

-    154 

3.90 

750 

33900 

-.313 

-       46. 

7 

6.70* 

800 

34  64  0 

.427 

53. 

1 

8.  04* 

850 

34765 

.616 

147 

4.  19 

900 

3455  1 

1.06 

236 

4.  49 

1000 

3451  0 

1.58 

402 

3.  93 

1250 

34825 

3.97 

772 

5.  14 

1500 

344  4  0 

4.69 

1102 

4.26 

average 

3445  1 

avera 

ge 

3.74 

deviation 

327 

deviation 

.47 

per   cent  .95 

*Values  not    used    fo: 


per    cent         12.6 
the   calculation    of   the   average 


Secondly,  as  the  frequency  gets  further  from  resonance  the 
cotangent  function  gets  larger  and  the  errors  also  get 
larger.  It  was  for  these  reasons  that  some  of  the  data  in 
Tables  VII-IX  were  deleted  for  the  calculation  of  the 
averages. 

This  completes  the  discussion  of  the  standing 
wave  tube  technique.  The  final  results  will  be  examined 
later  in  this  chapter  along  with  those  for  the  Chung-Blaser 
dual-microphone  technique.  One  additional  experimental 
technique    needs    addressing.       The      resonant    frequency    of   each 
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Helmhcltz  resonator  was  measured  independently  of  either 
method.  This  was  ione  using  a  signal  generator  anc  speaker 
to  excite  the  resonators.  Using  a  small  microphone  in  the 
resonators  to  detect  the  pressure  changes,  the  signal  ler.er- 
ator  was  tuned  to  the  frequency  which  resulted  in  the 
greatest      pressure      in     the        resonators.  These      resonant 

frequencies  were  549.5  Hz.  for  resonator  1,  398.5  Hz.  for 
resonator    2   and    795.0   Hz.    for    resonator    3. 

C.       CHUNG-BLASER    DUAL-MICROPHONE    TECHNIQUE 

The  dual-microphone  technique  developed  by  Chung  and 
Blaser  provides  a  quick  method  for  measuring  the  acoustic 
impedance  and  reflection  coefficient  in  the  laboratory. 
However,  as  will  be  shown  later,  the  technique  suffers  seme 
degredation  in  accuracy  compared  with  the  results  obtained 
by   the   5WT    method. 

1 .      Theory 

Chung  and  Elaser  have  shewn  [Ref.  1]  that  the 
complex  reflection  coefficients  at  the  two  microphone  loca- 
tions can  be  expressed  in  terms  of  the  transfer  functions  of 
the    incident  and   reflected    waves.    This    expression    is 


A®    =    -/££>-  (2    28) 

*fiZ     H/sn  {2'2) 


where  R;  (f)  and  Ra  (f)  are  the  complex  reflection  coeffi- 
cients evaluated  at  the  first  and  second  microphone  loca- 
tions, respectively.  H^*  ( f )  and  Hr(f)  are  the  acoustic 
transfer  functions  associated  with  the  incident  and 
reflected  wave  components,  respectively,  evaluated  at  the 
two    microphone    locations. 
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For  plane  wave  propagation  (no  mean  flew)  and 
neglecting   wall    losses,    H^-(f)    and   Hr{f)    can    be   expressed   as 

Hf.  (f)     =    9~jks  (2.25! 

and 

Hr  (f)     =    sJkS  (2.  30) 

where      "k"      is    the      waver.umber      and      "s"    is      the      microphone 
spacing. 

Application  of  equation  (2.28)  results  in  the 
reflection  coefficient  MR"  on  the  surface  of  a  sample 
material  not  at    a   microphone  location    being    expressed   as 

-  -  '-  -    =    -  —'-  -  -  (2     311 


where 

H»    =    e"^'  ,  (2.32) 


and 

H;    =    ejkl  (2.33) 

and    "1"      is   the      distance  from   the      first   microphone      to   the 
surface   of   the    sample   material.    Using    equations    (2.31-2.33), 

the    reflection    coefficient    at    the   sample    is 

B    ■    R,  <£UA/  (2.  34) 


where   [Ref.    1] 


R/  '   h&  -Hjn  (2*35) 


H,^(f)       is  the   acoustic  transfer   function   between    microphones 
1    and   2    and   H4(f)       and      Hr(f)       are   the    transfsr    functions   of 
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the    incident  ^nd  reflected    waves,    respectively,      between   the 
two    microphone    locations. 

Using   the    standard    relation    between    reflection    coef- 
ficient   end  complex    specific   acoustic   impedance 


(2.36) 


and  the  relationships  in  equations  (2.34)  and  (2.35),  the 
complex  specific  acoustic  impedance  can  be  determined  at  a 
distance   "1"   from   the   first    microphone  [ Ref .    1] 


(2.37) 


where 


(2.38) 


H    (f)    can   be  expressed   as 


(2.39) 


where  Ctl  and  Q  are  the  real  and  imaginary  parts  of  the 
cross-power  spectrum  of  the  signals  at  microphones  1  and  2. 
S ,,  (f)  is  the  auto-power  spectrum  of  the  signals  at 
microphone    1. 

Finally,      Seyfcert      and   Ross    show   (Ref.    15]      that   the 
phase   change  upon   reflection   is   given    by 


0(f)    =   tan-i 


(2.4  0) 
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where   Q^r     and  Ct-r    are  the      imaginary    and      real   parts      of   the 
cross-power  spectrum-of   the    incident    to    reflected    waves, 

2 •      Equipment  set-up 

Figure  2.7  shows  a  simplified  diagram  of  the 
Chung-Blaser  dual- microphone  technique  equipment  set-up.  At 
the  center  of  the  system  are  the  Hewlett-Packard  HF-35 
computer   and  HP-5420    Signal    Analyzer    System. 

a.  The    Hewlett-Packard   HP-85    Computer 

The  HP-85  computer  is  an  8  bit  microprocessor 
which  uses  the  BASIC  computer  language.  It  has  16K  bytes  of 
read/write  memory  of  which  14,579  bytes  are  available  for 
use  by  the  operator.  The  unit  used  in  the  study  had  its 
memory  expanded  to  32K  bytes  by  a  memory  module.  The  unit 
has  a  mass  storage  capability  in  the  form  of  a  magnetic  tape 
reader/recorder-  It  also  includes  a  127  millimeter  diagonal 
black  and  white  electromagnetic  deflection  CRT  and  a  32 
character    per      line    thermal    printer/plotter.  To    interface 

with    other      equipment,      an    I/O      ROM   and    interface      card   were 
added   to   provide   HP-IE    (IEEE-488)     interface    capability. 

b.  The    Hewlett-Packard    HP-5423    Signal   Analyzer 

The  Hewlett-Packard  HP-5420A  Signal  Analyzer 
provides  digital  time-domain  and  frequency-domain  analysis 
of  lew  frequency  analog  signals.  The  HP-542QA  system 
includes  the  HP-5447B  Digital  Filter  Module,  HP-544  10A 
Analog  to  Digital  Converter,  HP-5443A  Keyboard  and  Control 
Module  and  the  544  1 A  Display  Module.  The  equipment  used  in 
this  study  also  included  an  HP-10920A  HP-IB  interface  acces- 
sory. This  allowed  its  use  with  an  HP-IB  controller.  It 
allowed  the  transfer  of  data  and  measurement  set-up  informa- 
tion to  pass  to  and  from  the  controller  and  allowed  the 
reception   of  remote    control   commands.      Using    the    controller, 
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the      data      was      modified,       processed,         saved      and      recalled 
resulting  in  a    high    degree    of    flexibility. 

The  HP-542QA  Signal  Analyzer  perforrr.es  a  variety 
of  time- do  main  and  f reguency-domain  measurements.  Ml  meas- 
urements  are  done  using  digital  signal  processing.  This 
involves  filtering  and  sampling  the  input  waveforms.  The 
analog  input  is  first  filtered  by  a  lew-pass  filter  and  then 
passed  to  an  analog  to  digital  converter.  Here,  the  waveform 
is  sampled  which  results  in  the  conversion  of  a  voltage  (at  a 
certain  point  in  time)  into  a  numerical  representation. 
These  numbers  are  then  filtered  and  processed  digitally  to 
produce   the   various    measurements   preformed   by   the   analyzer. 

c.      The    Impedance   Measurement    Tube 

The  tube  used  in  this  part  cf  the  study  was  a 
piece  cf  3  inch  inside  diameter  PVC  pipe  12  inches  long  with 
a  3  inch  loudspeaker  mounted  at  oze  end.  Two  holes  were 
drilled  in  the  tube  and  fitted  with  7/8  inch  Tygon  tubing  to 
hold  the  two  microphones.  The  first  was  located  11. U2  cm. 
from  the  sample  end  (opposite  the  loudspeaker)  and  the 
second  was  located  4.59  cm.  from  the  first.  The  distances 
were    measured   to   the    center    of   the  microphones. 

The  microphones  used  were  General  Radio  1962  1/2 
inch  elect ret-condenser  microphones.  These  microphones  have 
a  sensitivity  level  cf  -40  dB  re  1  v/Pa  and  have  typically 
flat  response  curves  from  20  Hz.  out  to  1800  Hz.  which  cover 
the    frequency  range    of  the    study. 

3  •      Experimental   Procedures 

The  dual-microphone  technique  was  used  to  determine 
the  resonant  frequency  and  the  reactive  impedance  of  the 
same  three  Helmholtz  resonators  used  with  the  standing  wave 
tube  technique.  A  major  difference  between  the  two  techni- 
ques made  itself  known  almost  immediately.  The  time  to 
collect   and  reduce    the  data    was   vastly    different. 
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a.  Microphone   Calibration 

The    two    microphone/preamplif iar      channels   had   to 

be  calibrated  for  phase  and  amplitude  differences  prior  to 
use  for  the  actual  measurements.  Ih=  technique  used  for  this 
calibration  was  to  place  both  microphones  in  identical  sound 
fields  and  record  any  frequency-dependent  phase  and  ampli- 
tude differences.  This  calibration  data  was  then  used  by  the 
computer  to  condition  the  data  prior  to  the  calculation  of 
reactive  impedance,  phase  angle  and  complex  reflection  coef- 
ficient. In  this  manner,  the  effect  of  differences  between 
the  two  microphone  channels  on  the  data  was  minimized. 
Additionally,  this  method  also  allowed  for  imbalances  in  the 
A/D    converter   channels   of   the    signal    processor. 

The  identical  sound  fields  were  obtained  by 
mounting  the  two  microphones  at  the  same  radial  position  at 
the  sample  end  of  the  tube.  This  was  based  on  the  assumption 
that  the  sound  field  is  radially  symmetric  at  frequencies 
below  "rhe  sloshing  mode.  This  frequency  is  given  by 
[Ref.    18] 

4,=     o^c/d  (2-41) 

where    "  oc"    is   the   solution    to 

d£±2&  =o  (2.42) 

and  "d"  is  the  diameter  of  the  tube.  The  sloshing  mode 
corresponds  to  "«,"  which  is  equal  to  .5861.  This  results  in 
a  cut-off  frequency  of  2650  Hz.  The  calibration  program 
used    is    given  in   Appendix  A. 

b.  Validation  of  Tube   Performance 

After  calibration  of  the  microphone  channels, 
the  micrcphcnes  were  placed  in  the  normal  positions  in  the 
wall    of   the  tube   and   the      proper   operation   of   the    system   was 
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verified.  A  4.03  inch  long  piece  of  ?VC  pipe  with  the  sam  = 
diameter,  provided  with  a  rigid  reflective  end,  was  used  as 
a  sample  impedance.  The  expected  results  were  a  reflection 
coefficient   of    1    and    a   known   phase  change   equal   to 

Afm*gtmS#£  (2.»3, 

The  results  are  shown  in  Figures  2.3  and  2.9  Figure  2.8 
shows  that  the  reflection  coefficient  is  approximately  1  for 
frequencies  below  2.5  kHz.  except  at  1.5  kHz.,  were  the 
phase  goes  through  0  degrees  yielding  an  in determinancy 
which  is  discussed  in  section  e.  which  follows.  Above  2.5 
kHz.  non-plane  wave  modes  degrade  the  results.  The  measured 
slope  of  the  phase  shift  versus  frequency  shown  in  Figure  X 
is  -53  degrees  +  1  degree.  Using  equation  (2.41),  this  slope 
produces  a  sound  speed  of  356  +  1.34  m/sec,  which  is  cor.sxs- 
tant  but  significantly  higher  than  the  average  of  the  four 
values  obtained  by  the  S  WT  technique  listed  in  Tables 
II,VII-IX    (c   =    345.12   +    .38    m/sec). 

c.  Data   Collection 

Each  of  the  three  Helmholtz  resonators  was 
placed  at  the  sample  end  and  the  data  collected.  The  data 
was  tabulated  in  Tables  X-XII.  The  collection  and  recording 
of  the  data  was  done  completely  by  the  computer  and  the 
program    is    given    in    Appendix   B. 

d.  Data   Reduction 

After  the  raw  data  was  collected  the  first  step 
was  tc  determine  the  resonant  frequency  of  the  Helmholtz 
resonator.  The  data  collected  allows  this  to  be  done  in  two 
different  ways.  Beth  were  used.  First,  the  specific  reactive 
impedance   "J.s  //>c"      was    plotted   versus      the    frequency      and   a 
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TABLE    X 

■  ■                                   i 

Dual- 

Microphone 

Measured    Data 

Resonat 

or   #1    Run    1 

r« 

=    587.3 

+ 

2.2    Hz. 

ZZSZ- 

IsfrZ 

£& 

a  se    Angl 

a 

360 

♦   Phaser    A  no. 

Hz. 

Pa-s/m 

deg. 

deg. 

250 

-5.  87 

0 

360 

350 

-2.  83 

-    23.9 

336.  1 

400 

-1  .  91 

-    41.7 

318.3 

450 

-1.  23 

-    65.6 

294.  4 

500 

-.671* 

-105.2 

254.8 

550 

-.256* 

-147.6 

212.  4 

600 

.  1  14* 

166.2 

650 

.395* 

135.6 

700 

.6  76 

110.7 

750 

.919 

93.9 

1000 

1  .86 

55.2 

1250 

2.66 

43.  2 

1  500 

3.03 

36.  5 

Run    2 

f. 

=    590.6 

£ 

.2    Hz. 

Preg. 

h&Z 

£h 

a se    Ar.ql 

a 

360 

♦  Phase   Ang. 

Hz. 

Pa-s/m 

deg. 

deg. 

250 

2.  13 

19.2 

379.  2 

300 

-1.58 

-       8.4 

351.  6 

350 

-2.  14 

-    27.4 

332.6 

400 

-1  .  64 

-    52.  2 

307.  8 

450 

-1.07 

-    78.6 

281.  4 

500 

-.591 

-113.5 

246.  5 

550 

-.2  24 

-152.  6 

207.4 

600 

.0  76 

170.8 

650 

.  3  50 

140.3 

700 

.6  29 

114.  6 

750 

.8  90 

96.0 

800 

1.  19 

79.9 

900 

1.  52 

66.5 

1000 

2.22 

47.  8 

1250 

2.  96 

37.  1 

1500 

3.  17 

34.9 

_       _   j 
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linear   regression  was  used  done   using  the  two  points  or. 


either   side   o: 


w  h  e  r 


a  —  .a     t  h 


:  he   origin    ( 4      points    tot 3.1) 
tivs    impedance    was      equal  tc   zero.      This      gives   the    resonant 

frequency    sines    the    :?ac:-7e    impedance    of    a    resonator    aqu. 

zero    at  resonance.       The   graphs      and   resonant    frequencies    for 
each    resonator    are    shewn   in    Figures    2.10-2.14 

The  second  method  for  determining  the  resonant 
frequency  involved  the  use  of  the  phase  angle  between  the 
incident  and  reflected  waves.  Whan  the  phase  angle  was 
plotted  versus  frequency  (Figures  2.15-2.19)  resonance  was 
obtained  from  the  graph  at  the  poiat  where  the  phase  angle 
was  180  degress.  Again  a  linear  regression  was  used  with  the 
2  points  on  either  side  of  the  180  degree  phase  angle. 
Figures  2.15-2.19  also  show  the  fraquency  determined.  The 
resonant  frequency  used  in  subsequent  calculations  was  the 
average  of  the  two  frequencies  determined  from  the  two 
methods.  The  reproducibility  of  the  results  can  be  estimated 
from  the  ratio  of  the  standard  deviation  to  the  average  for 
the  four  (two)  separate  measurements  of  the  resonant 
frequencies  for  resonators  1  and  2,  which  were  0.5%  and 
2.0%,    respectively. 

Once  the  resonant  frequency  was  determined,  it 
was    usad   alcng    with    equation    (2.13) 


u  =  f  -   f./f 


(2.13) 


to  obtain  "X^".  Using  this  along  with  the  reactive  impedance 
ratio  "Xs//>c"r  the  experimental  value  for  the  ratio 


\M0 


(2.42) 
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■    —    —  ' 

TABLE    XI 

I 
j 

1 

Dual-M: 

.crop  hone 

Measured    Data  - 

Resonator   #2    Run    1           j 

i 

fa 

=    435.8    + 

3.8    Hz. 

££3S- 

£s#>S 

£h 

age   Angl^ 

360 

+    Pnass    Aricr. 

HZ. 

Pa-s/m 

deg. 

deg. 

200 

.215 

7.5 

367.5 

250 

-  4  .  2'4 

-       1.  6 

358.4 

300 

-.  151 

-       8.3 

351.7 

350 

-.436 

-    47.3 

312.7 

400 

-.206 

-121.7 

238.3 

U50 

.100 

164.  6 

500 

.3  86 

133.0 

550 

.6  67 

109.7 

600 

.  893 

94.8 

650 

1.07 

85.0 

700 

1.  27 

76.0 

750 

1.  40 

70.4 

800 

1.  60 

63.8 

900 

1  .  84 

56.7 

1000 

2.  12 

50.  1 

1250 

2.57 

42.2 

1500 

2.83 

38.  1 

Run    2 

■ 

■e 
— • 

=    421.7    + 

.  3    Hz. 

Frsjj. 

Zs&Z 

Pb 

ass    Anal 3 

360 

+  Phass   Ana. 

Hz. 

Pa-s/m 

deg. 

deg. 

200 

.  163 

1.6 

361.6 

250 

-2.  03 

-    41.7 

318.  3 

300 

-.9  67 

-    87.4 

272.6 

350 

-.4  34 

-132.5 

227.  5 

400 

-.0  82 

-170.5 

189.  5 

450 

.  171 

160.0 

500 

.  412 

134.4 

550 

.  5  93 

116.9 

600 

.813 

100.7 

650 

.9  93 

89.8 

700 

1.  27 

75.7 

750 

1.50 

67.  1 

800 

1.  78 

58.4 

900 

1  .  91 

54.8 

1000 

2.65 

41.2 

12  50 

3.35 

33.2 

1500 
i-_.  ._  

3.38 

33.0 

i 
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TABLE    XII 

\ 

Du 

al-fticrophcae 

Measured   Dat 

a   - 

Resonator    S3 

f,  =    94  7.  1    * 

1.0 

Hz. 

ILSSL- 

X$^2 

Phase    Anale 

360 

+   Phase    Ana. 

Hz. 

Pa-s/m 

deg. 

deg. 

400 

.719 

22.  1 

382.  1 

450 

.  1  74 

5.6 

365.6 

500 

-.0  39 

.9 

359.  1 

550 

-.5  97 

-    15.4 

344.  6 

600 

-1.01 

-    25.4 

334.6 

650 

-1  .  06 

-    34.6 

32  5.4 

700 

-.  921 

-    56.4 

30  3.  6 

750 

-  .  8  34 

-    69.4 

290.6 

800 

-.6  25 

-    94.6 

265.  a 

850 

-.3  80 

-126.9 

233.  1 

900 

-.201 

-150.7 

209.  3 

950 

.0  48 

173.5 

1000 

.  196 

155.2 

1050 

.3  35 

135.5 

1  100 

.571 

118.  1 

1200 

1  .05 

86.5 

1300 

1.34 

73.0 

1400 

1  .57 

64.7 

1500 
i  . 

1.  69 

61.2 

_    _      _        j 

was      calculated    for      comparison    with      the    theoretical      value 
obtained   from   equation    (2.23) 


**>  '      "c'AJ 


(2.23) 


Tables    XIII-XV       show    the   results    for       "K^p"    for    each      of   the 
resonators   tested. 

e.      Error   Analysis 

The  Chung-Biaser  dual-microphone  technique  had 
some  inherent  errors  built  into  it  by  virtue  of  the  equip- 
ment   construction  and   mathematical      techniques.       First,      the 
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TABLE    XIII 

! 

! 

! 
■ 

Dua i- A icr o p ho ne 

Reduced    Data    - 

Resonator 

#1    Run    1 

Frequency 

u 

ZzcP 

Hz. 

s-i 

X       10-3 

250 

-1142 

5.  14* 

350 

-    644. 2 

4.  39 

400 

-    470.0 

4.06 

450 

-    323.3 

3.80 

500 

-    196.0 

3.42 

550 

-       82.7 

3.  10* 

600 

20.0 

5.69* 

650 

1  14.6 

3.45 

700 

202.9 

3.33 

750 

286.0 

3.  21 

1000 

652.0 

2.  85 

1250 

971.6 

2.74 

1500 

1268 

2.39 

average 

3.  52 

deviation 

.79 

per    cent 

22.4 

Run    2 

Frequency 

*L 

£e*p 

Hz. 

s-i 

x     10-3 

250 

-1146 

1.  86* 

300 

-    863. 1 

1.  83* 

350 

-    646.9 

3.  31 

400 

-    472.3 

3.  47 

450 

-    325.4 

3.  29 

500 

-    197.9 

2.97 

550 

-       84.4 

2.65* 

600 

18.5 

4.  12* 

650 

113.2 

3.  09 

700 

201.  5 

3.  12 

750 

284.8 

3.  13 

800 

363.8 

3.27 

900 

512.  3 

2.  97 

1000 

651.  1 

3.  41 

1200 

909.2 

3.  26 

1400 

1151 

2.75 

average 

3.  17 

deviation 

.21 

per   cent 

6.  62 

♦Values   not    used 


for  the  calculation   of   tl 

ie   average. 

! 
_          j 
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TABLE    XIV  j 

Dual-Microphone    Beducefl   Data   -   Resonator    #2    Run    1 

Fr sgusnc^ 
Hz. 

200 

250 

30  0 

350 

400 

450 

500 

550 

600 

650 

700 

750 

800 

900 
1000 
1200 
1400 


IX 

Scip 

s-i 

X      10-3 

749.6 

.  29* 

509.7 

8.  32* 

33  3.  1 

.45* 

192.6 

2.  26 

74.8 

2.75 

28.0 

3.  58* 

120.2 

3.  21 

204.7 

3.  26 

283.  5 

3.  15 

357.  8 

2.99 

428.7 

2.  96 

496.8 

2.  82 

562.6 

2.  84 

689.0 

2.67 

810.  1 

2.62 

1042 

2.47 

1264 

2.  24 

aver 

age 

2.79 

ie viat 

ion 

.  33 

per    c 

ent 

11.8 

*X 

£«*/o 

s-i 

x    10-3 

689.6 

.  24* 

461.7 

4.  40* 

293.  1 

3.  37 

158.3 

2.74 

44.  8 

1.  83* 

54.6 

3.  13 

144.  2 

2.  86 

226.5 

2.  62 

303.5 

2.69 

376.3 

2.64 

445.  8 

2.85 

512.8 

2.  93 

577.6 

3.08 

702.3 

2.72 

822.  1 

3.  22 

1052 

3.  18 

1273 

2.  66 

average 

2.  91 

ieviat 

ion 

.  25 

Run    2 
Frequency 
Hz. 

200 

250 

300 

350 

400 

450 

500 

550 

600 

650 

700 

750 

800 

900 
1000 
1200 
1400 


per    cent         3. 59 
♦Values   not    used    in   ^he    calculation    of   the  average. 
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1 

TABLE    XV 

,  _  .   —   . .    _  1 

! 

Dual- 

•Hicroph 

cue 

R> 

aduced   Da1 

ta    -   Resona 

tor    *3 

Freauency_ 

*i 

^e*P 

Hz. 

s-i 

X      10-3 

aoo 

-  1844 

.  39* 

450 

-  1545 

50C 

-1295 

.03* 

550 

-  1082 

.06* 

600 

-    895.9 

1.  13 

650 

-    730.9 

1.  45 

700 

-    582.2 

1.58 

750 

-    446.8 

1.  87 

800 

-    322.0 

1.  94 

850 

-     206.0 

1.  84 

900 

-       97.3 

2.07 

950 

5.2 

9.  24* 

1000 

102.4 

1.  91 

1050 

195.2 

1.  97 

1100 

284.0 

2.  01 

1200 

452.0 

2.32 

1300 

609.  6 

2.  20 

1400 

758.9 

2.07 

1500 

901.6 

average 

deviation 

per    cent 

1.  87 
1.91 

.23 
12.0 

*   Values 

not 

use 

d     in 

the   calc 

illation   of 

the    average. 
..    _     _ 

construction  of  the  tube  limits  the  range  of  freguencies 
over  which  good  data  could  be  obtained.  The  diameter  of  the 
tube  determined  the  upper  limit  of  measurable  freguency. 
Since  the  tube  was  valid  only  for  a  plane  wave,  the  frequen- 
cies measured  had  tc  be  below  the  cutoff  freguency  for  the 
tube.  For  the  tube  used  in  this  study  the  cutoff  freguency 
was  approximately  2600  Hz.  The  construction  of  the  tube  also 
limits  the  effectiveness  of  measurements  at  the  lower 
freguencies.  The  critical  parameter  in  this  case  is  the 
spacing    between    the    microphones.       For      a    given   separation   as 
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the  wavelength  qets  longer  (frequency  lower) ,  the  differ- 
ences in  the  phase- and  amplitudes  at  the  two  microphone 
locations  b-caae  small  and  the  noise  dominated.  The  low 
frequency   limit    for    the   tube    was    about    200    riz. 

Finally,  the  mathematical  technique  resulted  in 
an  indeterminacy  in  the  reflection  coefficient.  This 
occurred  when  the  phase  angle  was  equal  to  zero.  When  this 
occurred   equation    (2.35)    reduced    to 

R    =   0/0 

which  is  indet  er  minant.  The  result  was  that  data  obtained 
for  the  frequencies  about  5  0  Hz  on  either  side  of  the  point 
where   the    phase    angle   equalled    zero    was    greatly    in    error. 

D.       COMPARISON    OF    THE    TWO   TECHNIQUES 

The  two  techniques  were  consistant  with  each  ether  and 
with  the  expected  theoretical  values.  First,  the  techniques 
will  be  locked  at  fcr  comparison  with  the  the  theoretical 
results   and   then    for    comparison    with    each   other. 

1 •      Standing   Wave  Tube    Technique 

Figures  2.20-2.22  show  the  comparison  of  the  experi- 
mental and  theoretical  values  for  "Xs//>c".  The  theoretical 
values  were  obtained  using  "Ky^"  (eqn.  2.23)  and  the  calcu- 
lated values  for  "X^"  (Tables  VII-IX) .  These  were  tabulated 
in  Tables  VII-IX  As  can  be  seen  in  Figures  2.20-2.22  the 
values  compare  favorably  with  a  divergence  noted  at  the 
higher  frequencies  away  from  resonance.  This  divergence  was 
explained   in  the   error  analysis   in   section   B.3.c. 

Table  XVI  shews  the  resonant  frequency  and  both 
constants  for  each  resonator,  along  with  the  resonant 
frequency  determined  by  the  independent  method  discussed  in 
section    B.3.c.         As    can   be    seen,      "&cXp"    ^.^d   "K/A"    were    within 


45 


the    calculated    error    for   each.  The   average   per    cent    error 

between  the  two  values  for  "K  "  was  found  to  be  '4,13;.  for 
the  standing  wave  tube  technique.  The  average  per  cent  error 
for   the    resonant   frequency    was    found    to    be    -2-23%. 

2  •      Chung- 3 laser    Dual-M  icrophone    Technigji  e 

Figures  2.23-2.27  show  the  comparison  of  the  experi- 
mental and  theoretical  values  for  ,tXs//°cn,  The  theoretical 
values  were  obtained  in  the  same  manner  as  in  the  SWT  tech- 
nique. The  graphs  show  a  good  agreement  except  at  the  lower 
and  higher  frequencies  away  from  resonance  for  reasons 
previously    noted   in    section    C.3.e. 

Table  XVI  shows  the  resonance  frequency  and  both 
values  obtained  for  the  dual-microphone  technique.  The 
values  determined  for  "K^,,"  and  "K^  "  do  not  in  all  cases 
overlap  within  the  areas  of  uncertainty  as  was  the  case  with 
the  SWT  technique.  As  a  result,  the  average  per  cent  error 
((K0<F-  K^)/K^)  x  100  was  found  to  be  13.4%.  The  average  per 
cent  error  between  the  independently  measured  resonant 
frequency  and  that  measured  by  the  dual-microphone  method 
was    found   to  be    +11.4??. 

3.      Standing   Wave  Tube    vs.    Dual- Micro£hone 

The  greatest  discrepancy  between  the  two  methods  was 
the  measured  resonant  frequencies.  The  resonant  frequency 
measured  by  the  dual-microphone  technique  was  always  higher 
than  the  resonant  frequencies  measured  by  the  other  two 
methods.  The  average  error  of  +11.4%  was  significantly 
larger  than  the  error  for  the  standing  wave  tube  method  of 
-2.23??.  Additionally,  the  determination  of  thfi  constants. 
for  the  dual-microphcne  method  produced  larger  discrepancies 
between  the  theoretical  and  experimental  values.  Finally, 
the  calculated  sound  speed  for  the  dual-microphone  method 
was       3.27c    higher      than      the    SWT      method      which      was    in      good 


46 


_ ._    . 

■ 

TABLE    XVI 

Cciisparisor 

i    of    Results 

Freq. 

K+K 

K£4P 

iktzzJU,  x  100 

Hz. 

sec. 
x    1  0- 

■3             X 

sec. 

10-3 

% 

RES   JM 

549.5   £    1.5 

3Z. 

SWT 

54  2.7 

7.35 

+    .  18 

7.57 

j£ 

.39 

+  2.99 

Dual   run    1 

58  9.  5 

3.6  7 

t    .  10 

3.52 

+• 

.79 

+  4.  09 

Dual    run    2 

59  0.  5 

3.66 

±    .10 

3.17 

± 

.21 

-13.4 

£ES    #2 

398.5   +    3.5 

Hz. 

SWT 

38  7.  4 

7.37 

±    -15 

6.94 

+ 

.41 

-5.83 

Dual   run    1 

43  2.  0 

3.45 

±    .33 

2.79 

+ 

.33 

-19.1 

Dual   run   2 

422.  0 

3.6  8 

j.    .07 

2.91 

i 

.25 

-20.9 

RES    #3 

795.0   ±    4.0 

Hz. 

SWT 

77  3.0 

3.6  1 

t    .09 

3.74 

_/■ 

.47 

+  3.60 

Dual   run    1 

97  4.  4 

1.42 

±    .04 

1  .91 

* 

.23 

+  34.5 

Hi    CENT    ER10R   II    i 
Resonator   #1    -   549. 

iESON&NI 

I    FRE£U 

ENCY 

5 

Hz. 

SWT         -      542 

.  7    Hz  . 

-  1. 

24% 

Dual      -      590 

.  3    Hz. 

+  7. 

4  2% 

Resonator  #2 

2   3  9  8. 

5 

HZ. 

SWT         -       38  7 

.4    Hz. 

-2. 

79% 

Dual      -      428 

.8    Hz. 

+  7. 

60% 

Resonator   #3 

-    795, 

0 

Hz. 

SWT        -      773 

.  0   Hz. 

-2. 

77% 

Dual      -       97  4 

.  4    Hz. 

+  22.6% 

• 

. j 
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agreement      with    the      expected      sound      speed,      based      on      the 
ambient    temperature. 
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Figure  2.  1    Helmholtz  Resonator  Physical  Set-Op. 
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Figure  2.2   Theoretical  Helmholtz  Resonator  Arrangement 
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III.     CONCLUSIONS 

Throughout  this  study  both  techniques  displayed  advan- 
tages and  disadvantages.  The  criteria  for  saying  one  is 
"better"  than  the  other  depends  upon  the  purpose  for  which 
the  information  is  needed.  No  attempt,  will  be  made  to  decide 
which  technique  is  "better".  The  strengths  and  weaknesses 
will  be  pointed  out  and  the  potential  users  of  this  informa- 
tion can  decide  which  technique  nore  closely  fits  -heir 
needs. 

It  was  clear  from  the  results  in  Table  XVI  that  the 
standing  wave  tube  technique  was  the  more  accurate  when 
compared  to  the  expected  theoretical  values  for  the 
Helmholtz  "standards".  In  addition,  the  SWT  technique  had  a 
higher  average  relative  precision  than  the  dual-microphone 
technique.  The  average  "for  the  SWT  was  4.2%  while  the 
average  relative  precision  for  the  dual-microphone  technique 
was  18.4%.  We  believe  the  reason  for  the  factor  of 
fcur-and-one-haif  better  results  for  the  SWT  technique  was 
due  to  one  major  characteristic  of  the  two  techniques.  The 
standing  wave  established  in  the  two  tubes  were  theoreti- 
cally the  same.  However,  the  locations  where  the  appropriate 
measurements  were  obtained  were  different  for  each  tech- 
nique. For  the  SWT  technique  the  measurements  were  always 
taken  at  the  optimum  locations,  ia.  a  maximum  or  a  minimum 
pressure  point.  Therefore,  it  was  possible  to  obtain  the 
best  data  available.  For  the  dual-microphone  technique  the 
measurements  were  taken  at  locations  on  the  standing  wave 
predetermined  by  the  microphone  locations.  These  locations 
were  only  optimized  when  the  microphone  separation  was  equal 
to  an  odd  multiple  of  a  quarter-wavelength  of  the  frequen- 
cies     measured.         Therefore,        the      data      obtained      was      not 
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necessarily  the  best  available  and  tha  result  was  a  degreda- 
tion  in  accuracy  and  relative  precision  for  the  iaal- 
ir.i cr o p hon  —  technique.  This  does  not  provide  an  obvious 
explanation  for  the  error  in  measured  sound  speed  or  signif- 
icant error  in  measured  Helmholtz  frequency  obtained  by  the 
dual-microphone    method. 

It  would  appear  that  the  SWT  method  was  the  more  desir- 
able based  on  accuracy  and  relative  precision.  This  may  not 
be  the  case.  When  numerous  samplss  were  involved  or  the 
results  desired  over  a  large  frequency  range,  the  SWT  method 
became  tedious  and  cumbersome.  It  involved  the  manual  taking 
and  reduction  of  data.  The  dual-microphone  technique  was  a 
great  deal  faster  and  easier  to  use.  The  time  to  obtain  the 
results  for  one  sample  over  a  frequency  range  of  200-1503 
Hz.  was  approximately  20  minutes  for  the  dual-microphcne 
technique    and   6    hours   for  the    SWT   technique. 

Thus  when  choosing  which  technique  is  more  desirable  for 
a  given  application,  it  must  be  decided  whether  speed  or 
accuracy  is  more  important.  The  degradation  in  the  accuracy 
and  relative  precision  for  the  dual-microphone  technique 
might  be  considered  minor  in  those  applications  where  speed 
is    important. 
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While  the  experiments  conducted  in  this  study  were 
straightforward,  several  changes  could  be  made  to  improve 
both  techniques.  The  biggest  change  would  be  to  adapt  the 
SWT  technique  for  computer  control.  This  could  be  accom- 
plished using  a  servo-mechanism  with  a  worm-drive  on  the 
microphone  car.  Ths  computer  could  then  automatically  set 
the  frequency,  determine  the  maximum  and  minimum  values  and 
their  locations  and  compute  the  necessary  information.  This 
would    greatly   enhance   the   SWT    technique. 

One  improvement  or  change  for  the  dual-microphone  tech- 
nique that  needs  investigation  is  the  use  of  a  pure  tone  as 
the  signal  instead  of  random  noise.  It  would  be  possible  to 
step  through  the  frequencies  desired.  The  use  of  a  pure  tone 
would  allow  the  use  of  a  narrow  band  filter  and  could  result 
in  a  reduction  in  the  noise  interference  and  would  include  a 
useful  redundancy  as  the  wavelength  and  thus  phase  differed 
between    the  two    microphones. 

Although  here  we  only  dealt  with  reactive  impedance,  it 
would  also  be  possible  to  extract  an  absolute  value  of  the 
resistive  component  of  the  Helmholtz  resonator  with  only  one 
additional  measurement.  If  the  quality  factor  of  the  reso- 
nator is  measured  and  the  absolute  reactive  impedance  is 
known,  the  absolute  resistive  impedance  is  fixed.  This  value 
can    then   be   compared    with   the    impedance    measurement. 

Finally,  several  changes  could  be  made  to  adapt  the 
dual-microphone  technique  for  use  in  a  water  filled  tube. 
First,  the  tube  should  be  constructed  with  numerous  hydro- 
phone locations  so  that  the  separation  between  hydrophones 
could  be  varied  somewhat  to  optimize  the  data  for  a  given 
frequency  range.      Additionally,    the   use    of    a   bubble   of    known 
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dimensions  could  be  used  to  determine  the  accuracy  ar.d  rela- 
tive precision  of  the  water  rube  system.  This  would  be  anal- 
ogous   to   the  Helmholtz   resonator    used   in   -his   study. 
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APPENDIX    A 
MICROPHONE    CALIBRATION    PROGRAM, 


1000 

1010 

102  0 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

option  base   1 

matrix    •»•    holds   HP5420    data    block    head=r 

matrix    'A'    holds   auto    and   cross    spectrum    data 

matrix    ' C*    is  calibration   vsctor 

address    of    HP-5420    is    704 
dim    W  (2,  16)  ,  A  (3,256)  ,C(2,256) 

output    704    ;"RS"   a    S=SPOLL    (704)     a)    wait    500    3    S#0    the: 
1060 

HP5420    test    setup: 

60    stable  averages 

channel    1    active 

transfer    function 

randon    signal 

1.6kHz    bandwidth 
output   704    ;"60r1    AV    CH    5    FR    2    SG     1.6    KV   BW    0    CH    0    AC 
10    RG" 

S  =  SPOLL    (704)     a   wait    500    a)    if    S  #)     then    1140 
if   S#0  then   1  1  40 
!    start    data    collection 

disp    "press   end  line    to   continue"    a   input    A$ 
output   704    ;"ST" 

S=SPOLL    (704)3    wait    5  00    5)    if    S#68    then    1190 
disp    "press   end  line    to   continue"    a)   input   A$ 
gosub    1330    a)    gosub    1600 

!    create    data    file'CALDAT'    if    not    present    on    current 
tape 

on   error    goto     1240   a)    create    "CALDAT" ,  1 ,  4096 
off   error    a)  assign*   1    to   "CALDAT" 
print*   1,1    ;    C  (,) 
!    close   file    to  record   data   on    tape 
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1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 


assigr.#    1    to    * 

disp    "calibration    data  store 


a  ii 


disp  "end  program" 

clear  7 

stop 

and 

!    collect,    data 

1    transfer    auto   spectrum    data    to    HP85 

output   704    ;"1    TC    1    FH    1    AU" 

S=SPOLL    (704)     2   wait    500    5)    if    S  #4    then    1360 

output   704    ;"501SA" 

S=SPOLL     (704)     a   wait    500    a    if    S#32    then    1380 

for   1=1    to    16 

!  enter  5420  header  information  into  matrix  w 

enter  704  ;  W  (1,1) 

next  I 

!  enter  data  in  matrix  A 

for  1  =  1  to  256  a  enter  704  ;  A  (1,1)  a  next  I 

!  transfer  cross  spectrum  data  to  H?85 

output  7  04  ;"RS" 

S  =  SFOLL  (704)  a  wait  500  a)  if  S#32  then  1470 

output  704  ;  "CR" 

S  =  SPOLL    (704)     a   wait    500    a    if    S  #0    then    1490 

output   704    ;"501SA" 

S=SPOLL     (704)     a   wait    500    a    if    S#32   then    1510 

for   1=1    to    16    a  ent=r    704    ;    »(2,I)     a   next    I 

for    1=1    to    256 

!    enter   real   cross   spectrum    data   to   matrix    'A' 

enter   704     ;    A  (2,1) 

!  enter  imaginary  cross  spectrum  data  to  matrix  'A' 

enter  704  ;  A  (3,1) 

next  I 

return 

!  calculate  calibration  vector 

print  "calculating" 
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1620  for  1=1  to  256 

1630  C  (1,I)*A  (2,1)  /A(1,I) 

1640  C  (2,1)  =A  (3,1)  /A(1,I) 

1 b 5 0  next  I 

1660  return 


82 


APPENDIX    3 
DOAL-MICBOPHONE    COMPUTER    PROGRAM 


10 
20 
30 
40 
50 
60 
70 
80 
90 

100 
110 
120 
130 

140 

150 

160 
170 
180 
190 
200 
210 
220 
230 
240 

250 


address    of   HP5420   si    704,    plotter    is    705 
list   of    variables 

mic   1    autospectrum 

cross   spectrum,    real    part 
cress   spectrum,    imaginary    part 
calibration   vector,    real    par- 
calibration    vector,    imaginary    part 
real   part  of   acoustical   transfer    function 
imaginary   part   of    acoustical   transfer 


a  v  "  r ) 
A{2,) 
A(3,) 
C(1,) 
C(2,) 
RH 
IH 

function 
A1  () 
PHI  () 

R() 
M() 


reflection    coefficient 
phase   change 
impedance,    real   part 
impedance,    imaginary    part 
option   base    1     a)   dim   A  (3  ,  256)  ,  A  1  (256)  ,  C  (2,  256)  ,  H  (256)  , 
R  (256)  ,PHI  (256) 

!    all   following  air   constants   and   calculations   are   in   SI 
units 

t  emperature 

specific   heat   ratio 

air    pressure 

gas   constant 

speed   of   sound   in   air 

air   density 


T3=295.15 

G3=1. 4 

P3=101325 

R3=287.0537 

C  =  SQR  (G3*R3*T3) 

RO  =  P3/(R3*T3) 


EOC=RO*C 

assign*  1  to  "CALDAT"  a)  read*  1,1  ;  C(,) 

disp    "enter   number   of    materials    to    be   measured" ;  2)  in  put 

18 

for  19=1  to  18  a)  print  "we  will  now  measure  material  #" 

;I9 
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260  disp  "if   sample    is    ready,    press    end   line"    3    input    AS 

270  cutout    704     ;  "RS11    S    S=SPOLL    (704)     S    wait    500    3    if    S#0 

ther.   270 
28  0  output    7  04     ;"30,1    AV    1    CH    5    FR    2    5G    1.6    KV    BW    0    CH    0     AC 

AC     10    RG" 
290  S=SPOLL     (704)    a    wait    500    3    if    5#0    then    290 

300  output    7  04     ;"ST    MR" 

310  S=SPOLL     (704)    5)    wait    500    9    if    S#68    then    310 

320  disp  "end   daxa   collection" 

330  output    704     ;"1    F1    1     IC    1     AU" 

340  S=SPOLL     (704)    3    wait    500    a)   if    S#4    then    340 

350  J=1    a)   gosub    530 

360  output    704     ;"RS"    a>    S=SPOLL    (704)     a)    wait    500    d    if    S#0 

then  360 
370  output    704     ;"CR" 

380  S  =  SPOLL     (704)    a)    wait    500    a)    if    S#0    then    380 

390  gosub    530 

400  output    704     ;  "RS"    5)    S=SPOLL    (704)     a)    waiz    500    5)    if    S*0 

then   400 
410  disp  "off   to    calculate!    this    will   take   several 

minutes" 
420  gosub    600 

430  disp   "press    end    line   to    store    impedance   data    on   tape" 

440  input   A$ 

450       !    the   following  statements    save    data    on    tape 
460      create  "IMPD AT" ,256, 2 4 
470       assign*    1    to    "IMPDAT"    a)    for    1=1    to    256    a)   F=(I-1)*D1    a) 

print#1,I;  F,  R  (I)    a)   next    I    d>   assign*    1    to   "IMPDAT2" 
480      disp   "end    data    collection   for   set";l9 
490      next    19 
500      disp    "end    program" 
510       end 

520       !    start    data    ccllection 
530       output    704    ;"501SA" 
540       S  =  SPOLL     (704)     a)    wait    500    5)    if    S#32    then    540 
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550   for  1  =  1  to  16  a)  enter  704  ;  H  o»  next  I 

560   if  J=2  than  580 

570      for   T  =  1   to   256    2)   enter  704    ;    MJ,I)    Snext    I    3>   print 

"auto-spectral   data  transferred"   3    J=J+1    2)    return 
580       for    1  =  1    to    256    3   enter   704    ;    A  (2, 1)     a)    enter    704    ;    A  (3,  I) 

dnext   I   a)    print   "cross-spectral    data   transf erred"areturn 
590       !    calculate   impedance 
600      disp    "calculating" 
610       for   1=1    to    256 

620  C1=C(1,I)*C  (1,1)  +C(2,I)*C(2,I) 

630  TA=  (A(2,I)  *C  (1,I)+A  (3,1)  *C(2,I)  )  /C1 

64  0  A  (3,1)  =  (A  (3, I)  *C  (1,1) -A  (2, 1)  *C(2,I)  ) /C1 

650  A(2,I)=TA 

660      next    I 
670       !    L   is  the   distance   from    the   sample    to    the    center   of 

the    farthest    mic. 
680       !    S   is   the  spacing   between    the    microphone   centers 
690       !    D1    is   the   frequency    spacing   of   the    5420    samples 
700       S=.0459 
710       L=. 1142    3    D1=6  .25 
720      !    calculate  transfer    function,    impedance,    reflection 

coefficient   and   phase    change 
730      for   1=1    to    256 
740  F=(I-1)*D1 

750  K=2*PI*F/C 

760  RH  =  A(2,I)/A  (1,1)     3    IH=  A  (3  ,1)  /A  ( 1  , 1) 

770  AA=SIN(K*L)     5)    BB=SI  N  (K*  (L-S)  )     a)    CC  =  COS  (K*L)     a)    DD=C0S 

(K*(L-S)  ) 
780  IN1=   BB-AA*RH   3    RN1=AA*IH 

790  RD1=RH*CC-DD    a)    ID1  =  CC*IH 

800  M2  =  RD1*RE1  +  ID1*ID1 

810  R(I)  ,RZ=  (RN1*RD1+IN  1*ID1)  /M2    a)    M  (I)  ,  I  Z=  (IN  1*RD  1  -RN1  * 

ID1)/M2 
820       RD=(RZ+1)/2 
830       ID=IZ/2    a)    M=RD*RD+ID*  ID 


8U0       RRH0--(RD/M-1)     3    IRHO=-(ID/M) 

350       &1  (I)  =SQR  <RRHG*RRKO+IRHO*IRHO)      3    PHI  (I)  =- AIN2     (IRHOfRRHO) 

373      next    I 

380       !    plot   results 

890       plotter    is    705    a   lorg    5    3    F0=500    3    F1  =  1600 

900       Y1=-10   3    Y2=10    a)    S  1=1 

910      disp    "resistance"   3   1=1    3   L$  =  "  RESISTIVE"   3    BS  =  " 

IMPEDANCE" 
920       gosub    1070    3    if   Z1    then    1020 
930       Y 1=- 10   9    Y2=10    a   S1=2 
940      disp    "reactance  a   1=2    a    L$="  REACTIVE"   a    B$=" 

IMPEDANCE" 
950       gosub    1070    a    if   Z1    then    1020 
960       Y1=0    a   Y2-1.4    a   S1=.2 
970      disp    "r-flection"    a   1=3    3   L$="  REFLECTION"    a    5$=" 

COEFFICIENT" 
980      gosub    1070    a    if  Z1    "hen    1020 
990       Y1=-130    a    Y2  =  180    a    S1=*45 
1000       disp    "phase   angle"    a    1=4    a    L3="  PHASE"    3    B$=" 

ANGLE     (DEG)" 
1010      gosub    1070    a    if   not   Z1    rhen   goto    430 
1020      disp    "lowest    Y"    ;a   input    Y1 
1030       disp    "highest    Y"    ;a   input   Y2 
1040      disp    "increments   of"     ;a    input    S1 
1050      on    I    goto     910,940,970,1000 
1060      !    plotting   subroutines 

1070      disp    "change    paper   and   press   end    line";ainput    A3 
1080       Y0=  (Y2-Y1)  /7/25   a    Y3=  (Y2-Y  1)  /7*  1.  9    a    Y4= (Y2- Y 1 0/7*1 .6 

a    Y5=Y1-Y3    a    Y6=Y2+Y4 
1090      lorg    5   a    scale    - (F1/5 ) , 1 . 16*F1 , Y5 , Y6    a    xaxis    Y1,F0rQ, 

F1    a    yaxis    0,S1,Y1rY2 
1100       Y7=Y1-Y3/6    a    Y8=Y1-Y3/2    a    ldir    0 
1110      for   X=0   to   F1       step   F0    3    move    X,Y7    a   label    VAL$ 

(X/1000)    a    next    X 
1120       move    .4*F1,Y8    a   label    "FREQUENCY-KHZ" 
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1130 

1140 

1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 

1230 
1240 


for   Y=Y1    to  y: 
(Y)    S  next   Y 


ldir    PI/2    S 


steo   3  1    a"  move   -(F1/20),Y   5)    label 
=  (Y2-Y1)/3  +  T1    a)    move    -(F1/10),Y9     J 


Y.-.  L$ 


L$,B$    a)    F9=0 

for    L=1    to   256    dtep   3 

on    I    goto    1170,1130,1190,1200 

F2  =  R(L)     3    goto    1210 

F2=M(L)     a)   goto    1210 

F2  =  A1 (L)     d    goto    121  0 

F2  =  PHI  (L)/PI*180 

move   F9,F2    a)    F9=F9+3*D1    d    iitove    -  (F1/500)  ,  YO/2 

idraw    F1/250,0   a)   iiraw    0,-Y0    a)   idraw    -(F1/250),)    a) 

idraw   0 , Y0 
next    L 
disp    "1    for   redraw,    0    to   continue"  ;a)   input    Z1    a)   return 
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